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(2) 

We present a toy model with a Hamiltonian Hf,’ on a folded one-dimensional spin chain. The 
non-trivial ground states of Hf ’ are separated by a gap from the excited states. By analyzing the 
symmetries in the model, we find that the topological order is protected by a Z 2 global symmetry. 
However, by using perturbation series and excluding thermal effects, we show that the Z 2 symmetry 

is stable in comparison to a standard nearest-neighbor Ising model with a Hamiltonian Hi. We 
(2) 

find that Hf is a member of a family of Hamiltonians that are adiabatically connected to Hi. 
Furthermore, the generalizations of this class of Hamiltonians, their adiabatic connection to Hi, and 
the relation to quantum error-correcting codes are discussed. Finally, we show the correspondence 
between the two ground states of H^ and the unpaired Majorana modes, and provide numerical 
examples. 


Intro- A large effort has been made to construct ro¬ 
bust protocols for quantum computation by exploiting 
topological properties of many-body systems [El [| . Mod¬ 
els such as the toric and surface codes |i| U are pro¬ 
posed. A particular attention has been given to detect 
and employ exotic non-Abelian excitations in quantum 
computation H, with a special attention to the Majo¬ 
rana fermions |j O . In this order, a crucial step is to de¬ 
tect and characterize the topological orders. It has been 
shown by Levin and Wen |t| that the ground state of 
quantum many-body systems with non-trivial orders can 
be seen as a condensate of fluctuating string-like objects. 
In particular, quantum phases are commonly studied by 
the projective symmetry groups (PSG) tool [|l3| which 
has been used to identify symmetry protected topolog¬ 
ical orders (SPT) [|t| [l]|, and with focus to quantum 
computation and quantum error correction HU- 
Consider a one-dimensional spin-1/2 Ising model with 
the Hamiltonian 

JV-l 

ib = -■/£«,, (1) 

i= 1 

on a chain with N spins (J >0) where a\ is the mth 
component of the Pauli matrices acting on site i. Ig¬ 
noring thermal excitation for a moment, a longitudinal 
field perturbation U z = lifts the degeneracy of 

the ground states for any non-zero hi. However, in the 
absence of a longitudinal field, the degeneracy is topo¬ 
logically protected aginst a perturbation U x = 

In comparison to the toric code [^) where the system 
is defined on a L x L lattice and the protection against 
perturbation is of the order of L, we could think of the 
one-dimensional Ising chain as a 1 x IV two-dimensional 
lattice. The lattice has a width of only “one” lattice site. 
Therefore, any (longitudinal) single-spin perturbation al¬ 
ready reaches the size of the lattice. Thus, the Ising chain 
has a topological phase, but the phase can not be real¬ 
ized due to the short width of the lattice. Although in 
different words, this claim was originally expressed in a 


footnote in Ref. [p3[ . 

In this Letter we introduce an adiabatic transforma¬ 
tion that does not change the topological characteristic 
of the Ising Hamiltonian Hi, but effectively folds the 
spin chain to a 2 x | lattice. As a result, we obtain 
a Hamiltonian Hf, ’ and the topological protection of the 
one-dimensional Ising chain tremendously improves, un¬ 
expectedly. In order to do so, the steps bellow are fol¬ 
lowed. 

/ON 

First, a Hamiltonian Hj ,' on an open chain of spin-1/2 
( 2 ) 

is defined. Hf, is a sum of four-spin operators. The set 
of these operators form a group that is denoted by . 
In addition, we find two symmetry groups <5>J 2 ^ and 
that commute with all the elements of . We find that 
(tS/ 2 ' 1 (g> 1 5>2 2 ' ) )/'!?( 2 ) = Z 2 which indicates the non-trivial 
order of the model is protected as long as the global sym¬ 
metry Z 2 is preserved. Next, we show through a de¬ 
generate perturbation analysis that the global symmetry 
Z 2 is robust. Then, the extension to models with wider 

width H^ vldth ' > and with symmetry groups {s / Wldth '>y i s 

( 2 ) 

discussed. Finally, we show the Hamiltonian Hf, is adi¬ 
abatically connected to Hi in Eq. ([]]), 

= R(ir) Hi (it) , (2) 

where R(a) is the unitary transformation 

R(a ) = e^ v . (3) 

Here V is a sum of two-spin interactions and a is a scalar 
parameter. The relation between Majorana modes of the 
Kitaev toy model [^3) and the two ground states of H^ 
are explained as well. 

The Model .- Consider a one-dimensional spin-1/2 sys¬ 
tem of a length N = 4m + 2, with a positive and integer 
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FIG. 1. (a) The standard one-dimensional Ising model with 
nearest-neighbor interactions, (b) The reshaped Ising model 
into a ladder after the unitary transformation R in spin space. 
The diagonal dashed and dash-dot lines show the two-spin 
terms crfcrj used in the transformation, (c) The Hamiltonian 

H {,' resulting from the transformation has four-spin interac¬ 
tions, associated to plaquettes. The bulk has two different 
types of plaquettes, A and B. There are three operators of 
type A and one operator of type B. The two left and right pla¬ 
quettes on the boundaries are named R and L, respectively. 
There are two operators acting on each of these boundary 
plaquettes. 


number m. We define H^ as 




-jEK’+f+f 


i=1 


— J (Li + R\ + Z/2 + R 2 ) j 


m —1 


i =0 

(4) 


where we have introduced the following operators (the 
stabilizers {O s }): 


A< i 1] = <i<i+i^li+^li+3 ( 5a ) 

A f ] = « + W4i+2 0'4i+3 ( 5b ) 

A f ] = <<+l°4i+2<+3 ( 5c ) 

B J = ^4j+l f7 4j+3 CT 4j+4^+6 ( 5d ) 

Ll = ^lO'2 cr 3 cr AT-l ( 5e ) 

(5f) 

L 2 = of of of (5g) 

R 2 = 0'fcr^_ 2 cr^_ 1 cr^, (5h) 


with i = 1,..., m — 1 and j = 0 ,..., m — 1. 

All the terms on the r.h.s. of Eq. (|dj) commute with 
each other. As shown in Fig. ilc, one can identify 
and Bi as bulk plaquettes operators, while R\ 2 and L\ >2 
act as boundary plaquette operators. For N = 4 m + 2 
there are 3(m — 1) plaquette operators of type A, m pla¬ 
quette operators of type B , and two plaquette operators 
of types R and L each. Therfore, there are overall 4m +1 


plaquette operators - that is N — 1. Without proof, the 
set of the stabilizers defined in Eq. (||) generates a group 
that we denote by . 

Separately, consider the set of the operators {of of}, 
on pairs {(*, J)}- Unless otherwise mentioned, we use the 
notation (i,j) for the two sites i and j that are con¬ 
nected with a dashed or a dot-dash line in Fig. [j]b. 
We define the group generated by these operators as 
S[ A = (of of,..., &N_ 2 a Ni a 2 a N-i)- AH elements in 
5-| 2) commute with all elements in . 

Furthermore, consider the set of operators {of of} on 
pairs {{i, j)}. We define the group generated by these op¬ 
erators as <S 2 2 1 = (of erf,..., ^N- 2 a N^ <j 2 (t n-i)- All ele¬ 
ments in <S 2 2 1 commute with all elements in and all 
elements in Especially, we find (iS{ 2 ' ) ® = 

Z 2 (there are N/2 generators in iS 2 2 \ AT/2 generators 
in s [‘ 2> . and N — 1 generators in <5>( 2 f 2 N /2 N ~ 1 = 2). 
A comprehensive classification of SPT orders in one¬ 
dimensional spin systems is given in Ref. Jl^| . 

The physical implication of the above statements is 
the following. If we manage to have a fixed “gauge” 
6 q d |G) = +|G) for the ground state subspace |G) of 
Hj, , as long as the Z 2 symmetry is not broken, the 
ground state is doubly degenerate and we can have |G±) 
such that 7 |G+) = e la / 2 |G_) for any generator 7 of <S 2 2 ^ 
fujf . The phase e l “/ 2 is a global phase and it is equal to 

e i7r / 2 £ or 7 /f 2 ) ( ] e fi ne( ] j n Eq. (|j). 

To see this, consider the states |f) = I"!') x ® • • • ® jt)^ 
and | J) = IJ.)! ® • ■ ■ < 8 > | 4 -)jv where by definition erf |t)j = 
+ |t) i and erf |J.) i = — Next, for every element Sk £ 

S[ 2) we have s k |t) z = + |f) z (»*= |X) Z = + |i) z ), and 
( 2 ) 

thus the gauge is set to +1. Then, by applying the 
group elements of S*- 2 ) on each of the states |f) and 
we obtain 

|G +>=y35n |i+c m w 

and 

^-^np+ojro.. < ? > 

where the product is over the stabilizers O s . Notice that 
2 to = A — 1 i s the number of distinguishable plaquettes 
in Fig. [yc. Thus, 4 m is the number of “loops” that can 
be constructed on the folded chain using the plaquettes 
as the unit blocks. The basis states appear as conden¬ 
sates of string-like configurations |12j. Interestingly, for 
a generator 7 = of of of <S 2 ' 

7 |G+) = i\G-), ( 8 ) 

since o y | {) = +i {). The equations (§, © , and the 
property in Eq. (^) can be examined in the examples 
given at the end of this Letter. 
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In fact, the property in Eq. (||) is not a coincidence if 
one remembers that by a Jordan-Wigner transformation 
the Ising model maps to the impaired Majorana problem 
pdf . Thus, the generators in the group S^ act equiva- 
lently as a logical operation for the ground states of Hf . 

Stability of the Z 2 symmetry In a one-dimensional 
Ising spin on an open chain a non-zero longitudinal field 
opens a gap between |f) and ||). The gap stimulates 
topological excitations (a propagating domain wall) from 
a false vacum to a true vacum and eventually destroys 


the symmetry 121 


( 2 ) 

Similar to Hi, the Hamiltonian Hf has a discrete 
energy spectrum and at low temperature (fc^T -C J ) 
thermal excitations are energetically costly and can be 
considered forbidden. By a degenerate perturbation se¬ 
ries approach we see that a transverse field perturbation 
such as U x = hi a i ( or U y = JA hio\) has vanishing 
matrix elements (G\U l x \G) for all the powers l < N, in the 
ground states supspace (|G) = |G+)(G+| + |G_)(G_|). 
Therefore, we have a topological protection with re¬ 
spect to this transversal field. Especially, the only non¬ 
vanishing term is {G\X = nj=i of |G). 

Considering the bases |G+) ± *|G_), we find X = 
nti °f an d any generator 7 G as the bit-flip and 
phase-flip logical operations, respectively. 

In contrast to Hi, in the presence of two transver¬ 
sal fields ( H^ +U X + U y ), the first non-vanishing term 
(G\U x U y \G) oc (G|of<rjf|G) appears at the second order 
of the perturbation series and only on pairs (i, j). This 
means that, out of (!)) number of possible terms in the 
second order, only ^ of them are non-zero. Thus, by in¬ 
creasing the length, the second order non-vanishing terms 
are suppressed by a factor of This is opposite to 

the Hi + U z case where in the second order of perturba¬ 
tion there are (^) — (N — 1) non-vanishing terms (O(N)). 
This pattern continues in all the orders. The odd orders 
vanish. In the fourth order there is a suppression factor 
of O(^), etc. 

Thus, as long as the perturbation (noise) affects sin¬ 
gle spins (no correlated noise) we should expect that the 
interplay between multiplicity and energy cost in the per¬ 
turbation series (the statistical ground for a phase tran¬ 
sition) to be substantially suppressed in our model with 
Hf Hamiltonian. Therefore, we expect a stable global 
symmetry Z 2 . Whether the Z 2 symmetry will still be 
destroyed through other mechanisms is a question to be 
answered. For a ladder of the toric code this has been 
studied Q. 


Translational symmetry breaking and generalization- 
By a close look at Fig. and Fig. |l|c, we notice 
that the periodicity changes as we move from the Ising 
chain in Fig. 0a to HIfi' 1 in Fig. 0c. The unit cell in 
HIfi' 1 is the two-plaques AB and the periodicity goes as 
• • • ABAB ■ ■ ■ in the bulk, while in the Ising chain it 


goes as • • • ZZZZ ■ ■ ■. The emergence of non-trivial topo¬ 
logical orders in one dimensional “organic” polymers by 
breaking translational symmetry has an old history [ ph| . 

Also, notice that the width of the folded chain in Fig. 
0c contains two sites. It should be now clear to the reader 
why we chose N = 4 x m+2. There are 4 sites (2 x width) 
in each unit cell and the last 2 sites are added to keep 
the inversion symmetry and have (iSj 2 '*®^ 2 ' 1 )/^*' 2 ^ = ^ 2 - 
Although, it is not clear whether the inversion symmetry 
is necessary. 

(3) 

Similarly, we can construct a Hamiltonian Hf ; with 
a width of “three” sites and define stabilizers with “six” 
operators ^(|. In this case, the length of the chain is 
chosen to be N = 6 m with a positive integer m. That is, 
to — 1 unit cells in the bulk. There are five different A- 
type operators, defined on A plaquettes, and one H-type 
stabilizer (in a hexagonal shape), with two boundary sta¬ 
bilizers on each edge. Therefore, there are AT — 1 number 
of stabilizers. The non-trivial part is to show that there 
are two independent symmetry groups S{ ; and S)) ' that 
commute with the group of stabilizers of H^. That 

is, they satisfy (<S^ 3 - ) ® S^)/S^ = Z 2 . 

(3) 

Hf' is the first member of the quantum error- 
correcting codes that can be constructed by folding a 
line and has a code distance of 0(1); That is roughly 
half of the width (3) of the folded configuration. In prin¬ 
ciple, it should be possible to construct quantum error- 
correcting codes with longer code distance and larger sta¬ 
bilizers with Hamiltonian H^f ldth \ As we will see for the 
case of Hf bellow, there is an adiabatic connection be¬ 
tween this class of Hamiltonians and Hi in Eq. (0). 

( 2 ) 

The Hamiltonian Hf ‘ is adiabatically connected to 
the one-dimensional Ising Hamiltonian Hi — We begin 
with the standard nearest-neighbor Ising Hamiltonian in 
Eq. (0) with J > 0. Next, the unitary transformation R 
defined in Eq. (||) with 

v = 5>r°7 ( 9 ) 

(id) 

is used to map Hi to H^ according to Eq. (||). The 
result is 

H(a) = R(a) Hi R^ (a) 

= «* 2 (f) ff, + cos(|)sin(|) H, 

+ sin 2 (|) 4 2 >, (10) 

(2) 

where Hi and Hf ; involve three-body and four-body in¬ 
teraction terms, respectively. For a = i r, the contribution 
of H\ to the total Hamiltonian drops out and we obtain 
H(n) = H^\ 

(2) 

Since we obtain Hf ; by a unitary transformation from 
Hi, the energy spectrum must stay unchanged and the 
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( 2 ) 

ground state subspace of Hp.' must be two-fold degen¬ 
erate. However, if two gapped states are connected by 
a set of local unitary transformations they belong to the 
same phase ^j. This can be seen by applying R(n) on 
a state |f)- One obtains 

|^)=iZ(7r)|t) = n(lt)ilt) J --*|;) j |l) J -), (ID 

(*.j) 

which is a product state. As it can be checked, clearly 
a second degenerate state is not accessible by using Eq. 
(||). That is, for a generator 7 of the group , we 
have 7 |V’) = — \^) while N = 4m + 2 . This means | ip) = 

(|G+) — i |G_))- This is not surprising if one thinks of 
the Majorana counterpart of the problem where the two 
unpaired Majorana modes (|G+) and |G_)) are paired 
up and experimentally undetectable. 

Then, how can we observe the phase that corresponds 
to |G+) (or |G_))? One quick answer is to find a way 
to implement logical quantum gates. We are seeking an 
operation such that M (± ) | ip) = |G±). Since 7 is a logical 
phase-flip operation and X = n)=i a f is the logical bit- 
flip operation, we should be able to decompose M (“) and 
have 

M (_) = -^= [<yX - i 7 ]. (12) 

Notice that the logical operation M is unitary but 
non-local. Implementing if ever possible experi¬ 

mentally, would change the quantum phase from a prod¬ 
uct state in Eq. ( 0 ) to a Majorana mode |G_). This 
can be checked for a chain with N = 6 in the examples 
bellow. 

Notice that we can obtain 7 and X, corresponding 
to the logical operations of Hp, , by transforming the 
X = nL <7f and the order parameter corresponding 
to Hi under R. In general, to obtain the logical oper¬ 
ations corresponding to H < j, nldth ' > one needs to know the 
adiabatic transformation from Hi and the knowledge of 
the symmetry groups of the Hamiltonian becomes irrele¬ 
vant. 

Examples and numerical results .- Let us consider two 
examples. The first example is a chain of N = 6. Fol¬ 
lowing Eq. (Q), we have (see Fig. |b) 

H^ n = 6 = -J(L 1 +L 2 +B 0 + R 1 +R 2 ). (13) 

By exact diagonalization we obtain the two basis states 

1 g+, n=6 ) = -J=[i tttttt) - 1 umt) - 1 mm) 

-1 mm)]- 

|g_, jv =6 > = -J=[| mm) - 1 nun) -1 umt) 

-1 nun)] 

(14) 



FIG. 2. Chains with N = 4 (a) and N = 6 (b) spins. The solid 
lines represent the initial Ising nearest-neighbor interaction 
between the spins, Hi. The dashed and dotted lines indicate 
the interaction terms present in V (Eq. (|^)) which are used 
to generate the transformation R. 


The second example is a chain with N = 4. It does 
not exactly follow the prescription defined in Eq. (Q). 
However, we can define (see Fig. ||a) 


#Tiv=4 = - J + of ufufer! + °203 oj), 

(15) 

and the inversion symmetry in this case is still preserved. 
It has the degenerate basis states 


iG +iiV =4) = ^[inn)-um)] 

ig_, iv= 4 ) = ^ [i uu) +1 mt)] ■ (i6) 

Notice that by applying ofo-f £ S^ on |G^) (|CG^ 6 )) 

one obtains VG^) (i\G^)). 

We use the dependence of the splitting of the ground 
states, A (a) = E+(a) — E_(a), on the global external 
magnetic field as a criterion to numerically verify the 
enhanced protection in Ht■ The numerical calculation 
of the splitting for N = 4 at the points a = 0 (Hi) and 
a = tt (Ht) is shown in Figs. and ||b as a function of 
h z (longitudinal) and h y (transverse) external magnetic 
fields, respectively. One can see that the dependence of 
the splitting A z goes from linear for Hi to quadratic for 
Ht, indicating increased protection. Figure ^b shows the 
gap A y which behaves topologically protected for both 
Hi and Ht, as expected. 

Summary.- The standard one-dimensional Ising chain 
with a Hamiltonian Hi could theoretically be in a topo¬ 
logical phase if the global Z 2 symmetry were stable. The 
symmetry is not stable since the order parameter of the 

system is just a single spin. We showed how to adiabati- 

( 2 ) 

cally obtain a Hamiltonian Hj ,' by a local unitary trans¬ 
formation with only two-spin interactions from Hi. We 
showed that the protection against single-spin errors in 
the transformed Hamiltonian Hip 1 scales with the length 
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FIG. 3. The energy splitting between the two low-lying states 
at a = 0 (Hi, circles) and a = 1 (Ht, squares) as a function 
of external field for a N = 4 spin chain, (a) The splitting as 
a function of a longitudinal field h z is shown for the two Hi 
and Hj. ’. For any non-zero longitudinal magnetic field h z , a 
gap opens linearly for Hi. The situation is visibly different 
(quadratic) for H^p as discussed by perturbation analysis, 
(b) The splitting as function of a transverse field h y is the 
( 2 ) 

same for both Hi and Hj , ; . 

of the chain. We discussed a family of Hamiltonians that 
are adiabatically connected to Hi. 

Acknowledgments - The author thanks Eduardo Muc- 
ciolo for his advise and support and Alioscia Hamma for 
his critical comments. This work was supported in part 
by the National Science Foundation grant CCF-1117241. 


[1] S. D. Sarma, M. Freedman, and C. Nayak, Phys. Today 
59, 32 (2006). 

[2] A. Kitaev, Ann. Phys. 321 , 2 (2006). 

[3] E. Dennis, A. Kitaev, A. Landahl, and J. Preskill, J. 


Math. Phys. 43 , 4452 (2002). 

[4] S. B. Bravyi and A. Y. Kitaev, arXiv:quant-ph/9811052. 

[5] A. Yu. Kitaev, Ann. Phys. 303, 2 (2003). 

[6] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. 
D. Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

[7] J. Alicea, Phys. Rev. B 81, 125318, (2010). 

[8] L. Fu and C. L. Kane, Phys. Rev. Lett. 102, 216403 
(2009). 

[9] E. Ginossar and E. Grosfeld, Nat. Comm. 5, 4772 (2014). 

[10] A. Biihler, N. Lang, C. V. Kraus, G. Moller, S. D. Huber 
and H. P. Biichler, Nat. Comm. 5 , 4504 (2014). 

[11] C. V. Kraus et al, New J. Phys. 14 113036, (2013). 

[12] M. A. Levin and X. G. Wen, Phys. Rev. B 71 , 045110 
(2005) 

[13] X. G. Wen, Phys. Rev. B 65 , 165113 (2002). 

[14] F. Pollmann, E. Berg, A. M. Turner, and M. Oshikawa 
Phys. Rev. B 85, 075125 (2012). 

[15] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Phys. 
Rev. Lett. 59, 799 (1987). 

[16] X. Chen, Z.-C. Gu, X.-G. Wen, Phys. Rev. B 83, 035107 

( 2011 ). 

[17] X. Chen, Y. Lu, and A. Vishwanath, Nat. Comm. 5 , 3507 
(2014). 

[18] W. Son, L. Amico, R. Fazio, A. Hamma, S. Pascazio, and 
V. Vedral, Europhys. Lett. 95 , 50001 (2011). 

[19] A. Miyake, Phys. Rev. Lett. 105 , 040501 (2010). 

[20] D. V. Else, I. Schwarz, S. D. Bartlett, and A. C. Doherty, 
Phys. Rev. Lett. 108 , 240505 (2012). 

[21] A. Altland and B. D. Simons, Condensed Matter Field 
Theory (Cambridge University Press, 2008), Ch. 3. 

[22] Notice that C S 1 ' 2 3 ' 1 . Thus, only one generator is 
enough to represent the whole subgroup. 

[23] A. Y. Kitaev, Phys.-Usp. 44 , 131 (2001). 

[24] V. Karimipour, L. Memarzadeh, and P. Zarkeshian, 
Phys. Rev. A 87 , 032322 (2013). 

[25] W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. 
Lett. 42 , 1698 (1979). 

[26] The details will be published elsewhere. 

[27] X. Chen, Z. C. Gu, and X. G. Wen, Phys. Rev. B 82, 
155138 (2010). 





